Application of cross correlations between CMB and large scale structure to 
constraints on the primordial non-Gaussianity 
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The primordial non-Gaussianity of local type affects the clustering of dark matter halos, and the 
planned deep and wide photometric surveys are suitable for examining this class of non-Gaussianity. 
In our previous paper, we investigated the constraint from the cross correlation between CMB lensing 
potential and galaxy angular distribution on the primordial non-Gaussianity, without taking into 
account redshift slicing. To improve our previous analysis, in this paper, we add the galaxy lensing 
shear into our analysis and take into account redshift slicing to follow the redshift evolution of the 
clustering. By calculating 81 power spectra and using the Fisher matrix method, we find that the 
constraint on the primordial non-Gaussianity can be improved from A/nl ~ 5.4 to 5.1 by including 
the galaxy-galaxy lensing shear cross correlations expected from the Hyper Suprime-Cam survey 
(HSC), in comparison with the constraint without any cross correlations. Moreover, the constraint 
can go down to A/nl ~ 4.8 by including the galaxy-CMB lensing cross correlations from the ACTPol 
and Planck experiments. 
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I. INTRODUCTION 

Primordial non-Gaussianities have been intensively 
discussed, because any detection of them offers an im- 
portant window into the early Universe. Standard 
single-field slow-roll inflation models predict small non- 
Gaussianity [1-5], so that the detection of large non- 
Gaussianity may turn down the standard model and sug- 
gest physics beyond the standard model. The most pop- 
ular method to hunt for primordial non-Gaussianities is 
to measure higher order correlation functions of the cos- 
mic microwave background (CMB) anisotropics such as 
the the three-point correlation function (bispectrum) and 
the four-point one (trispcctrum). Non- vanishing signals 
of these higher order correlations may predict the pres- 
ence of primordial non-Gaussianities. [6-11]. 

Recently, many studies have revealed that the primor- 
dial non-Gaussianity of local type affects the large-scale 
structure (LSS) through the clustering of dark matter 
halos, and shown the modification of the halo mass func- 
tion [12-17] and the halo bias [18-24], both by numer- 
ical simulations and by analytic calculations. Among 
several types of primordial non-Gaussianities, the local- 
type non-Gaussianity induces strong scalc-dcpcndcnce 
of the halo bias [18-21]. This scale-dependent bias 
is also powerful tool to constrain the primordial non- 
Gaussianity from the observations of LSS independently 
of the method with CMB higher order correlation func- 
tions. 
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Using this scale-dependent feature in the halo bias, 
some measurements of the primordial non-Gaussianity 
have been done in Refs. [18, 20, 25]. It is expected 
that future wide and deep surveys, such as Subaru Hy- 
per Suprime-Cam (HSC) survey [26], Dark Energy Sur- 
vey (DES) [27], Large Synoptic Survey Telescope (LSST) 
[28], and so on, will put tighter constraints on the pri- 
mordial non-Gaussianity, which are expected to be com- 
parable to those from the future CMB observations or 
more (e.g., [29-32]). 

To get a tighter constraint on the primordial non- 
Gaussianity, it is important to distinguish the signature 
of the primordial non-Gaussianity from the other effects 
more precisely, especially from the linear bias in the lim- 
ited survey area. Thus, it is better for this purpose to 
combine unbiased observables which are sensitive to the 
matter power spectrum in the near Universe, such as 
weak gravitational lensing of CMB and galaxy lensing 
shear. 

Cross correlations with complementary probes are ex- 
pected to provide additional information on top of their 
respective auto correlations. The cross correlation be- 
tween CMB and LSS is known to give information 
about the integrated Sachs- Wolfe (ISW) effect of CMB 
temperature anisotropy, which generates the secondary 
anisotropics due to the time variation of the gravitational 
potential [33]. Moreover, CMB lensing also has a corre- 
lation with LSS since the gravitational lensing of CMB is 
induced by the gravitational potential produced by LSS, 
and the galaxy lensing shear as well [34]. 

Since the effects of the CMB lensing are imprinted on 
small scales, the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite has not detected the CMB lensing di- 
rectly. Recent ground based CMB experiments such as 
Atacama Cosmology Telescope (ACT), announced the 
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detection of the CMB lensing signal, and an ongoing 
CMB observation by Planck [35] or various ground based 
experiments are expected to detect this signal more pre- 
cisely. Therefore, we can be sure that the observation 
of CMB lensing will increase its accuracy and play an 
important role in observational cosmology in near future 
surveys. 

In this paper we consider prospects for constraining the 
primordial non-Gaussianity through the scale-dependent 
bias from the near future surveys taking into account 
all auto and cross correlations among CMB, galaxy dis- 
tribution and galaxy lensing shear. We shall show that 
galaxy-CMB lensing and galaxy-galaxy lensing cross cor- 
relations are particularly fruitful and allow us to improve 
the constraints on the primordial non-Gaussianity. 

This paper is organized as follows. In section II 
we briefly review the effects of the primordial non- 
Gaussianity for large-scale structure and show the mod- 
ifications to the mass function and bias of dark matter 
halos. In section III we summarize the observables em- 
ployed in our analysis and describe the angular power 
spectra of the cross correlations between CMB and LSS. 
In section IV we mention the survey design and pho- 
tometric redshift systematics employed in our analysis. 
In section V we explain the method and the setup of 
our analysis. In section VI we show the results particu- 
larly focusing on the constraints on the primordial non- 
Gaussianity. In section VII we discuss the effects of the 
photometric redshift systematics and the massive neu- 
trino on the constraints. Then we assess the contribution 
of the cross correlations for the constraints. Moreover, 
we compare the constraints from some different surveys. 
Finally, in section VIII we summarize our conclusions. 
Through this paper we assume a spatially flat Universe 
for simplicity. 



II. PRIMORDIAL NON-GAUSSIANITY 



In the presence of the primordial non-Gaussianity, the 
mass function of clustering halos is modified and we 
adopt a non-Gaussian correction factor of the halo mass 
function based on the Edgeworth expansion [12]: 
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where S3 is the skewness of the density field which is pro- 
portional to /nl, ctm = o'm{M, z) is the rms of the linear 
density field smoothed on mass scale M, ^ is defined as 
V = 5cI(Jm and 8^ ^ 1.68 is the critical linear overden- 
sity. dn/dM is the mass function in the non-Gaussian 
case and duQ/dM is the one in the Gaussian case. For 
the Gaussian one, we adopt a model of Warren et al [38]. 

Recent studies have shown that the local-type pri- 
mordial non-Gaussianity produces a scale-dependent en- 
hancement of the clustering of halo on large-scales. 



Pg(A/, z, k) = b^iM, z)P{k, z) 



[6G+A6(M,z,fc)]2p(fc,z), 
(3) 



where Pg{AI, z, k) is the galaxy power spectrum, P(fc, z) 
is the matter power spectrum, bQ{AI,z) is the bias in 
the Gaussian case and Ab{M, z, k) is the non-Gaussian 
correction of the halo bias described as [24] 
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Here, fim.o and Hq are the matter energy density nor- 
malized by the critical density and the Hubble parameter 
at present, D{z) is the linear growth rate normalized to 
the scale factor a in the matter-dominant era and T(fc) 
is the transfer function of linear matter density fluctua- 
tions. For the halo bias in the Gaussian case bG{M,z), 
we assumed the form presented by Sheth et al [39]. 



Deviations from Gaussian initial conditions are com- 
monly parameterized in terms of the dimensionless pa- 
rameter /nl and the primordial non-Gaussianity of the 
local form is described as [3, 36, 37] 

$ = + /nl(<^'-(0')), (1) 

where $ is the curvature perturbation and is a Gaussian 
random field. On the subhorizon scale, 4> is related to the 
Newtonian gravitational potential 4* as $ = — 

The existence of the primordial non-Gaussianity, 
/nl 7^ 0, indicates that the initial density field is pos- 
itively or negatively skewed. Furthermore, the fact 
that the non-Gaussianity affects the clustering of dark 
matter halos or galaxies allows us to constrain on the 
non-Gaussianity from large-scale structure surveys. In 
particular, the local-type non-Gaussianity described by 
Eq. (1) induces a scale-dependent enhancement of the 
halo/galaxy power spectrum. 



III. ANGULAR POWER SPECTRUM 

Here we bricfiy review the auto and cross correlations 
of various cosmological observables and their angular 
power spectra. In this paper we take into account CMB 
temperature (T) anisotropics, E-mode polarization (E) 
and CMB lensing potential (-0) for CMB observables and 
galaxy distribution (g) and weak lensing shear (7) for LSS 
observables. 



A. Integrated Sachs- Wolfe (ISW) effect 

The decay of the Newtonian potential due to the pres- 
ence of dark energy produces a differential gravitational 
redshift, and this effect is called the late-time integrated 
Sachs- Wolfe (ISW) effect. In a fiat Universe the presence 
of ISW effect is a signature of dark energy, and induces 
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a non-vanishing cross correlation between CMB temper- 
ature and large-scale structure measurements, such as 
galaxy distribution, weak lensing field and so on. 

The contribution of the ISW effect to the CMB tem- 
perature field 0(n) can be written as 



6isw(«-) = 



ATisw(ri) 



CMB 



Tc 



CMB Jo 



dx^(xn,x), 



(5) 

where Tcmb is the mean temperature of the CMB, h 
is the direction to the line of sight, x is the comoving 
distance and x* denotes the distance to the last scatter- 
ing surface. $ is the gravitational potential and a dot 
denotes a derivative with respect to the conformal time. 

The angular power spectrum of the cross correlation 
between CMB temperature through ISW effect and the 
other measurements X can be written as 



(6) 



where 



Sa^oi^o^^ .z-|-^|T(fc),.(fcx(.)), 

(7) 

P<i{k) (X k" ^ is the primordial power spectrum of $ 
as a function of the wave number k and rf is the tilt 
of the primordial power spectrum. The functions T(fc) 
and D(z) are the transfer function and the growth rate 
for linear matter density fluctuations, respectively, and 
ji{kx) is a spherical Bessel function. The kernel Af{k) 
is for the other measurements, namely, CMB lensing po- 
tential, galaxy distribution and weak lensing shear in this 
paper; X = '0,g,7, respectively. 

The noise spectra of CMB include detector noise and 
residual foreground contamination. Here, we assume the 
ideal condition that foreground contamination can be 
completely removable and include only Gaussian random 
detector noise of the form [40] 



T,P 



|(6'fwhmAt,p)' 
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where ^fwhm is the spatial resolution of the beam and 
At,p represents the sensitivity to the temperature and 
polarization per pixel, respectively. These values are 
given for each of the frequency bands of channels v and 
we show the values for two CMB experiments such as 
Planck [35] and ACT with new polarization sensitive re- 
ceiver (ACTPol) [41] in Table I. 



B. CMB lensing potential 

CMB photons are deflected due to the gravitational 
potential produced by the large-scale structure on the 
way propagating to us (e.g., Ref. [42]). The relationship 
between the lensed temperature anisotropy T{Ti) and the 
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/sky 


V 


^FWHM 


At 


Ap 






[GHz] 


[arcmin] 


[/iK /pixel 




Planck 


0.65 


100 


9.5' 


6.8 


10.9 






143 


7.1' 


6.0 


11.4 






217 


5.0' 


13.1 


26.7 


ACTPol 


0.05 


148 


1.4' 


3.6 


5.0 



TABLE I: The specifications for CMB experiments, /sky is the 
sky coverage, Sfwhm is the beam width at FWHM, At and 
Ap represent the sensitivity of each channel to the tempera- 
ture and polarization, respectively, and v means the channel 
frequency. 



unlcnscd one T(fi) is given by the deflection angle d(n) 
as Tin) — Tin -f d). The deflection is related to the line 
of sight projection of the gravitational potential \E'(n, z) 
as din) = \/tp(n), where 



dx 



X*X 



*(xn,x). 



(9) 



Here, ip{'ri) is the (effective) lensing potential. 

The angular power spectrum of the lensing potential 
can be written as 



PdkP^{k)[AJ{k)] 



(10) 



where 



X.-xnk)D{z{x)) . , 



X*X a(x) 

(11) 

The lensing potential can be reconstructed with 
quadratic statistics in the temperature and polarization 
data that are optimized to extract the lensing signal. 
To reconstruct the lensing potential, one needs to use 
the non-Gaussian nature imprinted into the lensed CMB 
statistics, and the noise of the lensing potential can be 
-1 estimated as the reconstruction error [43, 44]. In this pa- 
^ per, we estimate the noise spectrum of lensing potential 
Nf following the technique developed in [44] optimally 
(8) combining the temperature and polarization fields. 



C. Galaxy distribution 

The luminous sources such as galaxies must be the 
most obvious tracers of the large-scale structure in the 
linear regime, and the projected density contrast of the 
galaxies can be written as 

f°° nAz) 
5g,i{n)=l dzbcff{z,k)^-^d{x{z)n,z), (12) 

Jo 

where the subscript i represents the i-th redshift bin. 
S{xn, z) represents the matter density fiuctuation, and 
ni{z) and are the galaxy redshift distribution and the 
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total number of galaxies per steradian in the i-th red- 
shift bin. The function 6eff (z, k) is the weighted effective 
halo/galaxy bias defined as 



besiz,k) 



dM- 



dn 
dM 



dMb{M, z, k) 



dn 



Ma: 



dM' 
(13) 

where Mobs is the observable mass threshold, which is 
the minimum mass of the galaxy we can observe, and we 
take the value to be Mobs = 10^^-'^[h~^MQ]. The func- 
tion b{AI,z,k) and dn/dM are the halo mass function 
and the halo bias in the case of the non-Gaussian initial 
condition defined in Eqs. (2) and (3), respectively. We 
show the dependence of 6off(z, k) on the redshift and the 
wave number with different mass thresholds and /nl in 
Fig. 1. 

The angular power spectrum of the galaxy distribution 
between i-th and j-th redshift bins can be written as 



where 
Af(fc) = 



k^dkP^{k)Af{k)Af{k), (14) 



dzbM^,k)'^T{k)D{z)je{kx). (15) 



To estimate the signal-to-noise ratios and errors in pa- 
rameter determination for each survey, we need to de- 
scribe the noise contribution due to the finiteness in the 
number of sources associated with source samples. We 
can write the noise spectra from the shot noise as 



N: 



giSj 



(16) 



where n,; is the mean surface density of sources per stera- 
dian in the i-th redshift bin. 



D. Weak lensing shear 
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FIG. 1: The effective weighted bias defined in Eq. (13) as a 
function of redshift z (Top) and as a function of wave number 
k (Bottom). (Top) We plot the Gaussian case (/nl = 0) 
for the values log Mobs ~ 11.2,11.7,12.2, respectively, and 
the k dependence does not appear in this case. (Bottom) 
We plot the non-Gaussian case (/nl ~ 50, 100) at redshift 
z = 0.5, 1.0, 2.0 as indicated, and the value log J\fobs is fixed 
to be log Mobs = 11.7. 



The weak lensing shear (equivalently the convergence 
in the weak lensing limit) is a weighted integral of the 
density field of sources, which is directly related to the 
source galaxy redshift distribution. The average conver- 
gence of a light ray bundle from sources in the i-th red- 
shift bin is written as 



(17) 



where Wi (x) is the convergence weight function of kernel 
defined as 



W^ixiz)) 
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Here xi^) is the distance to the lens and x(^') is the one 
to the source. 



The angular power spectra of the weak lensing shear 
between i-th and j-th redshift bins can be written as 



HI 3 



k^dkP^{k)Aj^{k)A]'ik) 



(19) 



Af{k) = / dzW,{z)T{k)D{z)Mkx). (20) 



where 



The measurement of the shear from galaxy images has 
some uncertainties, one of which mainly comes from the 
intrinsic shape ellipticities of galaxies. The galaxy shapes 
can be treated statistically through a shape noise contri- 
bution and the noise spectra can be written as 



an 



(21) 
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where is the intrinsic galaxy shear. In this paper, we empirically derived value [45] for all surveys, 
assume the intrinsic galaxy shear to be cr^ = 0.3 as the 

I 



survey 




rig [arcmin ^] 


survey area 






DES 


0.8 


10 


5,000 deg^ 


(/sky 


~ 0.12) 


HSC 


1.0 


30 


2,000 deg^ 


(/sky 


~ 0.05) 


LSST 


1.2 


50 


20,000 deg^ 


(/sky 


~ 0.50) 



TABLE II: Survey parameters adopted in this paper for DES, HSC and LSST. The parameter jZm denotes the mean redshift. 
The parameter fig denotes the mean number density of source galaxies at all redshifts, where the source galajcies are divided 
into five redshift bins in all surveys. The ranges of redshift are summarized in Table III. 
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< 2 < 











TABLE III: The redshift ranges of the i-th bin in the case of dividing into five redshift bins. is the mean redshift for each 
survey shown in Table II. We assume five redshift bins and adopt these redshift ranges for all surveys. Zmax is the limitation 
of observations and we assume Zmax = 4.0 for all surveys in this paper. 



E. Cross correlation angular power spectrum 



The effect of the primordial non-Gaussianity stands 
out the most in the galaxy-galaxy auto correlation (gg) 
through the bias parameter. As many previous works 
show, almost all the constraints on the primordial non- 
Gaussianity come from the galaxy-galaxy power spec- 
trum [18, 20, 25, 46, 47]. The cross correlation power 
spectra, e.g., CMB Temperature-galaxy cross correlation 
(Tg), however, provide extra information and improve 
the errors of the parameters. 

In our previous analysis [48], we paid a particular at- 
tention to the CMB lensing-galaxy cross correlation (TAg) , 
and estimated the contribution to the constraint on the 
primordial non-Gaussianity. Here, we improve our anal- 
ysis by adding the galaxy lousing shear, which may also 
correlated with galaxy distribution and CMB lensing, 
and estimate the errors of determining cosmological pa- 
rameters including any possible auto and cross correla- 
tions. 



Fig. 2 shows an example of the cross correlation angu- 
lar power spectra, for galaxy-galaxy lensing (g7), galaxy- 
CMB lensing (g-0) and CMB lensing-galaxy lensing ("07) ■ 
The boxes around the curve show the expected measure- 
ment errors for HSC-Zifce survey (for an illustrating pur- 
pose, we take the logarithmic binning, as ~ 0.23£). 
For the cross correlations with CMB lensing, we show the 
errors for two cases that one is for Planck and the other 
is Planck with ACTPol. 



IV. SURVEY DESIGN 

The errors in determining cosmological parameters 
from observations of LSS depend a great deal on the sur- 
vey design, such as the survey region, the mean redshift 
and photometric redshift errors. Here we summarize the 
survey design employed in our analysis. 



A. Redshift distribution of galaxies 

We assume the redshift distribution of galaxy samples 
with a function of the form: 



/3 z° 



(22) 



where is the mean number density of source galaxies 
at all redshifts, a, /3 and arc free parameters. We 
adopt a = 2.0, /3 = 1.5, and the parameter zg is related 
to the mean redshift Zm as 

ng r[(a + l)//3] 

and zp is determined in such a way that the mean redshift 
Zm fits the value in Table II for each survey. The nor- 
malization of the redshift distribution function is fixed 
by the total number density of galaxies: 

/•oo 

= / dzn(z). (24) 
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FIG. 2: The angular power spectra for galaxy-CMB lensing {gi^; Top), galaxy-galaxy lensing (gi7j; Middle), and CMB 
lensing-galaxy lensing (i/'Ti; Bottom). Here the subscripts i, j represent the redshift bin and we divided into 5 redshift bins 
for galaxy distribution (g) and cosmic shear (7). The angular power spectra shown here are i = j = 1,3,5 from left to right, 
respectively. The boxed bars are the errors of the power spectra assuming HSC, Planck and/or ACTPol, and we adopt the 
logarithmic binning for illustrative purpose. We also plot the angular power spectra in the non-Gaussian case, /nl ~ 50, 100 
and MobB[h~^ Mq] = 10^^'^, 10^^'^. For comparison, we plot the contribution of the photo-z error, ai''' = 0.03 for each bin. 



B. Photometric redshift systematics 

In the galaxy imaging survey whose observables arc 
galaxy distribution and galaxy shear, we measure a large 
number of galaxies so that the systematic uncertainties 
play a correspondingly important role with statistical er- 
rors. Although systematics may include many effects and 
these effects depend on the details of the observations 
[49, 50], we include only the effect of photometric red- 
shift errors. 

Imperfect calibration of photometric redshifts induces 
a residual scatter and bias with respect to the true red- 
shifts, and uncertainties in the redshifts distort the vol- 
ume element. Following the photometric redshift model 
as described in [51], the probability distribution of pho- 
tometric redshift Zph given the true redshift z, p{zpii\z) 
is modeled as a Gaussian distribution 



Pizph\z) 



1 



■ exp 



^ph ^bias) 
2^2 



(25) 



where crz{z) and Zbias(2:) are the scatter and bias, respec- 
tively, which are arbitrary function of redshift z. 

The best-estimated distribution for objects in the i-th 
photometric redshift bin with z^^ < z < z^^^\ ni(z), 
can be written as 



ni{z) 



-ph 



dzphniz)p{zph\z), 



= -n(z) [erf(xj+i) - erf(a;j)] , (26) 

where Xi is defined as 

X, = (zg -z + Zbias(^)) /V2a,{z). (27) 

In the same way as Eq. (24), the total number density of 
galaxies in the i-th bin becomes 



dzni{z). 



(28) 
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We assume a fiducial redshift scatter and bias following 

[50] as 

a,(z) = aW(l + z), (29) 
^biasW = 4iL(l + ^)' (30) 

where a^f' and z^jl^^ are defined for each redshift bin. We 

set the fiducial values as ci''' = 0.03, z\^}^^ ~ for each 
redshift bin and all surveys. 



CMB LSS Cross: [CMBxLSS] 



CP 


TT, EE, ^^j, TE, 


gg, 77, g7 


Tg, T-y, ?/)g, ii-y 


p 


3000 (8000) 


800 


800 


/sky 


0.65 (0.05) 


0.05 


0.05 



TABLE IV: Survey parameters for our analysis. We assume 
Planck (ACTPol) for CMB experiments or their combination. 
For the LSS survey, we assume the HSC survey and the sur- 
vey area of HSC fully overlaps with Planck or ACTPol exper- 
iments. 



FORECASTS 



• LSS + Planck + ACTPol: 



We estimate the parameter errors using the Fisher ma- 
trix formalism. We summarize the Fisher matrix for- 
malism taking into account the cross correlation between 
CMB and LSS, and we describe the setup for the CMB 
and LSS obscrvables. 



17 /Cross \ ^ 77I 

— /sky 2^ -^i 



■Cross I 
y I ACTPol 



/•ACTPol 

/sky 2^ -^i 



■CMB I 
ij I ACTPol 



A. Fisher matrix formalism 

Under the assumptions that the observables follow the 
Gaussian statistics, we can obtain the information on a 
set of parameters p from the Fisher matrix written as 
[52] 



//■LSS fCrossN T-Cr 
Usky /sky ) 



OSS 

«i I Planck 



Planck 
sky 



/sky^) X! 



CMB I 
iJ I Planck 



— ./sky 



<'max 

E 



(2^ + i; 



Tr [Ce,,C-^C,,,Ci^] , (31) 



where /sky is the sky coverage of the experiment, Ce is 
the signal plus noise covariance matrix, and Cg^i is the 
derivative of C with respect to the parameter pi. 

For CMB experiments, we consider two cases: that one 
is the constraint from Planck and the other is from Planck 
with ACTPol. It is well known that combining the re- 
sults of the satellite CMB experiment with the ground 
based one can improve the constraints, e.g., combining 
WMAP with ACT[53], ACBAR[54] or so on. For these 
two cases, we define the total Fisher matrix as below, 
respectively, under the assumption for the size of survey 

QvoQ Qc fACTPol ^ yLSS ^ fPlanck tTovo fACTPol j-LSS 

area as j^^^ < /^j^y < j^^^ . Mere j^^^ , /^j^y 

and /iky "'^'^ are the sky coverage of the ACTPol, the LSS 
and the Planck surveys, respectively, and we define the 
overlap region between CMB experiments and LSS sur- 
veys as /i^™"^ 

• LSS + Planck: 



where 



Cross 



/ j.Planck _ /-ACTPol \ x-CMB I 

[J sky / sky ) f-,ij I Planck ' 



(33) 



-/ I 



CMB 

,i3 



{21+1] 



(2^ + i; 



j'^Cross ^^Cross ^ ~ ^ ^Cross ^^Cross ^ ~ ■'■ j 

(34) 

Tr [QCf B(<^CMB)-lc.CMB(<^CMB)-l] ^ 

(35) 



and different superscripts for ^max and /sky represent dif- 
ferent observations or CMB experiments, e.g., i'^^ and 
^max^ are the maximum multipole we can take for the 
CMB experiment and the LSS survey, respectively. Here 
^min is the minimum multipole we can take and we use 
^min = 2 in all cases. We summarize these values in Table 
IV. 

The covariance matrix Ci and each of its matrix ele- 
ments is defined as 



Tp /-Cross \ ^ x-Cross i /-Planck \ ^ x" 

— Jsky 2-^ "T /sky /-^ < 



CMB 
,i3 



„in=2 



I / /-Planck 
+ Usky 



/iLSS 



/sr) E (32) 



JxY = Cj, + Nf. & 



XY, 



(36) 



where X,Y = 



T,E,'>p,g,,^, (for Cross) 



(for CMB) 



Cf^ is the angular power spectrum and Nf-^ is the noise 
spectrum for each auto or cross correlation. We defined 
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Nf^ and iVf ^ as Nj and N[ in Eq. (8), Nf^' and 

N]^'^^ in Eqs. (16) and (21), respectively, and Nf'^ is the 
noise spectrum of CMB Icnsing potential we estimate fol- 
lowing the technique developed in [44] . Then we assume 
that there is no correlation of noises between the differ- 
ent observables and different redshift bins. Moreover, 
we assume that i?-mode polarization [E) has correla- 
tion only with temperature anisotropics (T) because E- 
mode polarization can be generated through the Thom- 
son scattering dominantly on the last scattering surface. 

Incidentally the assumption that C^'^ = is not 
strictly correct [55] because the large-scale polarization 
i5-mode generated by scattering at reionization correlates 
with sources of CMB lensing potential over long distance. 
Moreover, if the galaxies are at high- redshift, the assump- 
tion that C^^' = is also not correct [56]. However the 
influence of these assumptions for our results, especially 
constraints on the primordial non-Gaussianity, can be 
small because sample galaxies in our analysis are at low 
redshift. 

Here we explain the detail of our Fisher matrices given 
in Eqs. (32) and (33). In Eq. (32) the first term is respon- 
sible for the cross correlation between LSS and CMB, the 
second one is for CMB observables at the higher multi- 
poles than the LSS survey can probe, and the last one 
is for CMB observables in the remaining survey area of 
CMB which does not overlap with the LSS survey. In 
Eq. (33) the first term is responsible for the cross corre- 
lation between LSS and CMB with ACTPol, the second 
one is for CMB observables with ACTPol at the higher 
multipoles than the LSS survey can probe, the third one 
is for the cross correlation between LSS and CMB with 
Planck in the overlapping region which does not over- 
lap with ACTPol, the forth one is for CMB observables 
with Planck in the region which does not overlap with 
the LSS survey, and the last one is for CMB observables 
with Planck at the higher multipoles than the LSS sur- 
vey can probe in the region which does not overlap with 
ACTPol. 

In our analysis, the total number of non-zero angu- 
lar power spectra is 81 for the full cross correlation case 
{T, E,ip,gi,ji), while it is 5 for the only CMB case 
iT,E,^). 



B. Setup 

Given the measurement design of the surveys, we can 
estimate the errors in determining the cosmological pa- 
rameters using the Fisher matrix formalism. The for- 
malism tells us how well the given surveys can measure 
the cosmological parameters around fiducial cosmological 
model. Therefore, the parameter forecasts by Fisher ma- 
trix formalism depend on the choice of the fiducial model 
and the number of free parameters. 

As for our fiducial cosmological model, we assume a 



flat ACDM model and we include the following 9 cosmo- 
logical parameters, which is based on the WMAP seven- 
year results [57]. The density parameters of cold dark 
matter, baryon and dark energy are Qch^, ^hh^ and JIa, 
respectively; dark energy equation of state parameter is 
w; the optical depth at reionization is r; and the primor- 
dial power spectrum parameters are the spectrum tilt rig 
and the amplitude of the primordial power spectrum A|j 
which is normalized at ko ~ 0.002 Mpc^^. We also in- 
clude the primordial non-Gaussianity of local type mod- 
eled by the non-linear parameter /nl- For specifying the 
galaxy bias, we include the minimum mass of halo host- 
ing the galaxies, for which we can observe Mo}^s[h^^ Mq] 
as a free parameter. The fiducial values are 

{lOOn^h^ n,h^, Qa, t, nS, w, /nl, logA/obJ 

{0.20, 0.1109, 0.72, 0.086, 0.96, 2.4xl0"^ -1, 0, 11.7}. 

Throughout this paper we assume a spatially flat Uni- 
verse, and the Hubble parameter is adjusted to keep our 
Universe flat when we vary Q\. 

To calculate the angular power spectra C; including 
non-linear effects on the angular power spectrum of the 
galaxy distribution, cosmic shear and lensing potential, 
we use the CAMB code [58] and the HALOFIT code [59]. 
We, however, want to remove the burden of the uncer- 
tainty of the non-linear calculation as much as possible, 
so our estimation preferably is based on the information 
from large-scale for LSS where the linear prediction is 
reliable. 

In our estimation, we include the information from 
temperature anisotropics (T), _E-mode polarization (E) 
and reconstructed lensing potential (V-') for CMB. The 
range of multipoles are 2 < £ < 3000 and the survey area 
is taken as /sky = 0.65 for Planck and /sky = 0.05 for 
ACTPol, respectively. On the other hand, for the galaxy 
survey (g) and galaxy weak lensing survey (7), the ranges 
of multipole are 2 < £ < 800 and the survey areas are 
/sky = 0.05 for both of them with HSC. Furthermore we 
take into account all the cross correlation we can con- 
sider and we assume the optimal condition that there is 
no correlation between different patches, so that the area 
having a correlation between CMB and LSS corresponds 
to the overlapped survey area. We summarize the values 
we use in the following calculation in TABLE IV. 



VI. RESULT 
A. Signal-to-Noise ratio 

To see the errors in the cross correlation power spectra 
more quantitatively, we estimate the signal-to-noise 
ratio (S/N) for each £-bin and the cumulative (S/N). 
We define the {S/N)'^ for a cross correlation between X 
and Y, following [60], as 
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FIG. 3: Signal-to-Noise ratio for some angular power spectra. We show the S/N for each bin (box) and stacked S/N (hne), 
respectively. Top panels are for galaxy-CMB lensing {gitp) angular power spectra, Center ones are for galaxy-galaxy lensing 
(gi7j), and Bottom ones are for CMB lensing- galaxy lensing based on the fiducial model of Fig. 2. 



S/N for each ^-bin: 



-j (£)^/3ky£(2^+l) 



{cr? 



• stacked S/N: 



(37) 



Jj] (<Cax)^/sky5](2^+l) 

1=2 



(Cf^ + NI^){CJ'' + iV/^) + (CfY)2 ■ 



(38) 



where ^max and are the maximum and mini- 

mum multipoles of the i-th multipole bin, respectively, 
Cj"^ and N]^^, NJ^ are the fiducial and noise spec- 



tra, respectively. The index X, Y represents {X, Y} G 

{r,£;,v,g,7}- 

Figure 3 shows an example of (S/N)"^ for each £ given 
by Eq. (37) and cumulative (S/N)^ by Eq. (38), for the 
galaxy-galaxy lensing (g7), galaxy-CMB lensing {gip) and 
CMB lensing-galaxy lensing (ip^) cross correlations. For 
the galaxy-CMB lensing and CMB lensing-galaxy lensing 
with HSC-Planck experiments, the amplitude of (S/N)'^ 
decreases at high-£ bins and the stacked {S/N)'^ satu- 
rates due to the noise of CMB lensing at small scales. 
On the other hand, we can expect the lager (S/N)'^ at 
high-£ bins for the case with HSC-ACTPol. The am- 
plitude of {S/N)'^ at low-£ bins where the signature of 
the primordial non-Gaussianity shines out is very low in 
all the cases. However, comparing galaxy-CMB lensing 
with galaxy-galaxy lensing signals, the amplitude of the 
former is larger than the latter in high-redshift bins and 
the difference becomes even greater for the with HSC- 
ACTPol. 



J 
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HSC + Planck + ACTPol 



CMB + g 

CMB + g + \|/ 

CMB + g + Y 

CMB + g + \\i + j 




1.04 -1 -0.96 0.72 0.73 0.74 0.75 11.68 11.7 11.72 



0.96 0.72 0.73 0.74 0.75 11.68 11.7 

£2. losM, 



FIG. 4: Projected lcr(68%) confidence constraints in some parameter spaces for the fv^O model. We sliow contours from 
CMB and galaxy distribution (black thick); from CMB, galaxy distribution and CMB lensing (red dashed); from CMB, galaxy 
distribution and galaxy lensing (green dot-dashed); and from all the observables (blue thin). All auto and cross correlations 
between these observables are taken into account for the constraints. Here we assume the HSC survey and Planck experiment 
(Left) or HSC and combination of Planck and ACTPol (Right). 





lOOOb/i^ 






r 








•w 


/nl 


log Mobs 


galaxy clustering only (HSC): Cf^ 


0.94 


0.026 


0.011 


0.11 


0.029 


0.13 


0.49 


0.053 


5.9 


0.60 


HSC + Planck 






















CMB only 


0.015 


0.0013 


0.19 


0.0050 


0.020 


0.0047 


0.042 


0.69 






CMB + Cf 


0.011 


0.00055 


0.010 


0.0044 


0.0034 


0.0027 


0.029 


0.042 


5.4 


0.027 


CMB + Cf + Cj" 


0.011 


0.00055 


0.010 


0.0044 


0.0034 


0.0027 


0.029 


0.042 


5.3 


0.027 




0.011 


0.00053 


0.010 


0.0043 


0.0033 


0.0027 


0.028 


0.041 


5.3 


0.027 


+ cr + cj' + ^ 


0.011 


0.00031 


0.0059 


0.0044 


0.0024 


0.0024 


0.027 


0.022 


5.1 


0.024 


+ c/'* + cT + cp + c;" 






















+ + cf ^ + cp 


0.011 


0.00031 


0.0059 


0.0043 


0.0024 


0.0023 


0.026 


0.022 


5.1 


0.024 


HSC + Planck + ACTPol 






















CMB only 


0.0076 


0.0010 


0.18 


0.0046 


0.018 


0.0038 


0.036 


0.66 






CMB + Cf 


0.0076 


0.00052 


0.010 


0.0041 


0.0033 


0.0025 


0.028 


0.042 


5.4 


0.026 


CMB + Cf + Cp 


0.0076 


0.00052 


0.010 


0.0041 


0.0033 


0.0025 


0.028 


0.042 


5.3 


0.026 


+ cp + cp + cp 


0.0076 


0.00050 


0.010 


0.0039 


0.0030 


0.0024 


0.026 


0.040 


5.0 


0.024 


+ cp + c;^ + cp 


0.0074 


0.00030 


0.0059 


0.0041 


0.0024 


0.0022 


0.025 


0.022 


5.1 


0.024 


+ cp + cj' + cp + C7» 






















+ + cp + cp 


0.0073 


0.00029 


0.0058 


0.0039 


0.0021 


0.0022 


0.024 


0.022 


4.8 


0.021 



TABLE V: Expected marginalized error (Icr) and these forecasts correspond to the results of Fig. 4. 



B. Constraints on the primordial non-Gaussianity 

We estimate the errors in the determination of cos- 
mological parameters for the HSC-Zifce survey using the 
Fisher matrix formalism, and show the main results, 



namely, the Icr confidence limit constraints in Fig. 4 and 
Table V, in which we consider 10 cosmological parame- 
ters: the fraction of the massive neutrino density parame- 
ter to the matter density parameter /iy(= 0.03) = r2j//57m 
in addition to the nine cosmological parameters defined 
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-5 



0.73 0.74 



FIG. 5: Projected lcr(68%) confidence constraints in some 
parameter spaces from Planck and HSC experiments includ- 
ing auto and cross correlations between CMB, CMB lensing, 
galaxy distribution and galaxy lensing. We show the con- 
tours with or without the photometric redshift scatter for the 
ACDM model, the contours without the photometric redshift 
scatter for the /i/ 7^ model and the contour with the photo- 
metric redshift scatter for the fv=0 model. 



in Sec. VB. 

The simplest case, which is from only galaxy clustering, 
puts the constraint on /nl as A/nl ^5.9, and the con- 
straint can be improved to A/nl ~ 5.4 by combining the 
CMB-only result though CMB has no information about 
primordial non-Gaussianity on the level of the 2-point 
function without cross correlations with galaxy cluster- 
ing. This is because the information from CMB helps 
to determine the galaxy bias by breaking degeneracies 
with other cosmological parameters, and then the degen- 
eracy between /nl and Mobs can be broken effectively. 
Hence the constraint on Mobs is also improved dramat- 
ically after combining CMB information, which can be 
seen in Table V, and these tendencies have been seen in 
our previous analysis [48]. 

We show four cases to see the improvement by adding 
the different observables in Fig. 4. In the first case we use 
CMB and galaxy distribution (CMB + g); in the second 
case we add CMB lensing to the first case (CMB -t- g + 
-0); in the third case we add shear to the first case (CMB 
+ g + 7)1 and in the last case we use all the information 
available. We take into account all the available power 
spectra. The selection of auto and cross correlations in 
Fig. 4 and Table V corresponds to each other in each 
case. We obtain the constraints on /nl as A/nl ^ 5.1 
with Planck and HSC. Moreover, we can improve the 
constraint by including ACTPol as A/nl 4.8, thanks 
to the galaxy-CMB lensing cross correlations. 



VII. DISCUSSION 
A. Selection of the fiducial model 

We assume imaging surveys for the LSS galaxy surveys, 
and we have to take into account the influence of the var- 
ious systematics, e.g.,, magnification effects due to gravi- 
tational lensing [61], photometric redshift(photo-2:) errors 
[51] and so on (see, for example, [49, 50]). Moreover, as 
mentioned in the previous section, the parameter fore- 
casts by Fisher matrix formalism depend on the choice 
of the fiducial model and free parameters. In Fig. 5, we 
summarize the effects of selecting fiducial models. 



1. Photometric redshift error 

For the constraint on /nl, any significant difference 
is not found between with and without the photometric 
redshift scatter systematics. This is because the redshift 
scatter changes the overall amplitude of the galaxy dis- 
tribution (g) and cosmic shear (7) as shown in Fig. 2. 
On the other hand, the effect of non-Gaussianity, /nl, 
emerges only at large angular scales, and it can be distin- 
guishable. However, the effects of the photometric red- 
shift error are significant for the parameters which are 
determined through the information about the redshift 
evolution such as JIa and w. Although the effect of the 
parameter Mobs is to change the bias parameter and is 
similar to the effect of the redshift scatter, a significant 
difference is not found in the figure. 



2. Massive neutrino 

Various experiments, such as ground based experi- 
ments of neutrino oscillation, predict the non-zero neu- 
trino mass [62, 63]. Therefore, we need to include the 
mass of neutrinos for a more accurate estimation of the 
cosmological parameters. There are, however, some de- 
batable subtle arguments about taking the mass of neu- 
trinos into account in models for the mass function and 
the halo bias. These models have been tested by N-body 
simulations without massive neutrinos and are not guar- 
anteed enough for cosmological models with massive neu- 
trinos. The analytic formula for the scale-dependent bias 
presented in Eq. (4) is also for models without massive 
neutrinos. 

For the mass function, it is found that the influence 
of massive neutrinos is very small [64]. Thus, we just 
assume that the halo mass function and halo bias models 
can be applied even if the massive neutrinos are included. 

Comparing the results oi = and ^ models 
in Fig. 5, a significant difference can be seen in the con- 
straint on Mobs- The existence of massive neutrinos af- 
fects the large-scale structure formation through its large 
velocity dispersion and alter the amplitude of the matter 
power spectrum. This effect is similar to the effect from 
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FIG. 6: Same as Fig. 4 but we highlight the contribution from the cross correlations to the constraints. The base constraints 
are from CMS{TT,EE,TE), galaxy distribution (gg), CMB lensing (iptp) and galaxy lensing (77) auto correlations (black thin 
line). The other contours show improvements by including cross correlation as indicated in the figure. Here we consider the 
HSC survey and Planck with and without ACTPol experiment (left and right panels, respectively). 
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X 10" 


w 


/nl 


log Afobs 
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CMB + cf + cf'^' + c;^ 


0.011 


0.00034 


0.0086 


0.0043 


0.0024 


0.0024 


0.026 


0.033 


5.4 


0.024 


+ cj' 


0.011 


0.00034 


0.0086 


0.0043 


0.0024 


0.0024 


0.026 


0.033 


5.3 


0.024 


+ ct^ 


0.011 


0.00034 


0.0086 


0.0043 


0.0024 


0.0024 


0.026 


0.033 


5.4 


0.024 


+ cp 


0.011 


0.00034 


0.0086 


0.0043 


0.0024 


0.0024 


0.026 


0.033 


5.4 


0.024 




0.011 


0.00031 


0.0059 


0.0043 


0.0024 


0.0023 


0.026 


0.022 


5.2 


0.024 


+ c1'> + c/^ + cp + Cf< 


0.011 


0.00031 


0.0059 


0.0043 


0.0024 


0.0023 


0.026 


0.022 


5.1 


0.024 


HSC + Planck + ACTPol 






















CMB + Cf + Cf'*' + Cf 


0.0074 


0.00033 


0.0086 


0.0039 


0.0022 


0.0022 


0.024 


0.032 


5.4 


0.021 




0.0074 


0.00033 


0.0086 


0.0039 


0.0022 


0.0022 


0.024 


0.032 


5.3 


0.021 


+ cf 


0.0073 


0.00032 


0.0085 


0.0039 


0.0022 


0.0022 


0.024 


0.032 


5.4 


0.021 


+ cf' 


0.0073 


0.00033 


0.0085 


0.0039 


0.0022 


0.0022 


0.024 


0.032 


4.8 


0.021 


+ 


0.0073 


0.00030 


0.0059 


0.0039 


0.0021 


0.0022 


0.024 


0.022 


5.1 


0.021 


+ cf " + cp + ct<^ + cr 


0.0073 


0.00029 


0.0058 


0.0039 


0.0021 


0.0022 


0.024 


0.022 


4.8 


0.021 



TABLE VI: Expected marginalized errors (Icr) and these forecasts correspond to the results of Fig. 6. 



Mobs and these two parameters are strongly degenerated 
with each other as shown in Fig. 4. Note that the de- 
grade of the constraint on Mobs in the case of the /zy 7^ 
model comes from the strong degeneracy between them 
and we have confirmed that there is little change in the 
constraint if we fix the value of fi, (and do not take it as 
a free parameter) . For the constraints on /nl itself, how- 
ever, the difference of the fiducial models whether = 



or not is not so significant. 

In summary, focusing on the constraint of /nl, the 
selection of the fiducial model makes little impact on the 
result. On the other hand, the constraints on some other 
parameters depend on the fiducial model, especially on 
the existence of neutrino mass and photometric redshift 
error. 



B. Contribution of cross correlations ing more efficiently than Planck. 

Now let us discuss the contribution from ACTPol ex- 
periment which will extract the information of CMB lens- 
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HSC + Planck HSC + Planck + ACTPol 







f. 
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./nl 


log Afobs 






w 


./nl 


log Afobs 


CMB + Cf + + CJ~< 


0.015 


0.0028 


0.054 


5.5 


0.033 


0.015 


0.0023 


0.053 


5.4 


0.029 


+ cj' 


0.015 


0.0028 


0.053 


5.4 


0.033 


0.014 


0.0023 


0.053 


5.3 


0.029 


+ ct' 


0.015 


0.0028 


0.053 


5.5 


0.033 


0.014 


0.0023 


0.053 


5.4 


0.029 


+ ct' 


0.015 


0.0028 


0.053 


5.4 


0.033 


0.014 


0.0023 


0.051 


4.8 


0.028 


+ Cf^ 


0.0095 


0.0028 


0.035 


5.2 


0.033 


0.0095 


0.0023 


0.034 


5.2 


0.029 


All 


0.0095 


0.0028 


0.034 


5.1 


0.033 


0.0092 


0.0023 


0.033 


4.8 


0.028 



TABLE VII: Expected marginalized errors (Icr) which correspond to the results of Fig. 8 with all systematics. This table should 
be compared with TABLE VI where the systematics are neglected. 



As shown Fig. 4 and Table V, we find that almost all 
the constraints on each plane are determined through 
galaxy lensing information except for Mobs and /nl. 
However, CMB lensing also gives important contribu- 
tions if the ACTPol experiment is included. Focusing on 
the constraints on the /NL-ATobs plane, both the CMB 
lensing and galaxy lensing signals improve the constrains 
slightly by including their cross correlations, though the 
most of the constraints are determined by the galaxy auto 
correlation. 

We show how the cross correlations improve the con- 
straints in Fig. 6 and Table VI. We find that most of 
the constraints are mainly determined through the auto 
correlations on each plane but we can still gain a little 
benefit from the cross correlations. By including all cross 
correlations among CMB, galaxy distribution and galaxy 
lensing shear, we can expect 5.6% improvement for the 
constraint on /nl with HSC and Planck, and 11.1% im- 
provement with HSC, Planck and ACTPol. In particular, 
the impact of the cross correlation between galaxy and 
galaxy lensing (g7) is significant, especially in the w-Vtp^ 
plane. This result is consistent with that obtained by 
[31, 61] and g7 will play a very important role in the 
constraint on dark energy parameters. 

Finally, to see the impact of cross correlations for the 
constraint on the primordial non-Gaussianity, let us focus 
on the constraints on the /NL-Mabs plane in Fig. 6. Wc 
find that the cross correlation between galaxy and CMB 
lensing (gV') will improve the constraint on /nl when the 
Planck and ACTPol experiments are combined. However 
we have little benefit from the cross correlations if we use 
only Planck for the CMB experiment. 



C. Dependence on mass threshold A^obs 

The parameter A/obs reflects the mass threshold wc can 
observe as galaxy and relates to the galaxy bias through 
Eq. (13). We can interpret from Fig. 1 that the larger 
value of Mobs predicts the larger bias. We show the con- 
straints for the different fiducial value of Mobs in Fig- 7. 
Only in this figure, we plot A log Mobs instead of log Afobs 
for the purpose of illustrating the different fiducial mod- 
els on the same planes. 



The difference which comes from the different fiducial 
values of Mobs appears especially in the constraints on 
./nl and Mobs, while the constraints on the other param- 
eters are not altered so much. This is because the con- 
straints on /nl and A'/obs come mainly from the galaxy 
distribution (g), while the constraints of the others come 
mainly from the other observables. As for the constraints 
on /nl, the larger value of Mobs can lead to a tighter 
constraint on /nl- The larger value of Afobs means pick- 
ing out the higher mass objects. Because high mass ob- 
jects exhibit a large bias, the non-Gaussian correction to 
the bias becomes large. Furthermore, the correction to 
the mass function also becomes large for high mass ob- 
jects. For these reasons, the effect of the non-Gaussianity 
through the effective bias &cff becomes clearer when we 
choose the larger values of Mobs for the fiducial model. 
Thus we can constrain /nl more tightly for larger the 
value of Mobs- 



D. Effects of photometric redshift systematics 

Here, we investigate effects of photometric redshift sys- 
tematics on the parameter constraints following examples 
in [50]. 

In our analysis, we assume the fiducial redshift scatter 
and bias as Eqs. (29) and (30), respectively. To charac- 
terize the uncertainty in the scatter and bias, we param- 
eterize ffi'^ and 4iL = {!' 2, 3, 4, 5}) in 5 bins, include 
these 10 parameters in the Fisher matrix, apply a prior of 
0.01 on each, and finally marginalize over these param- 
eters. We show the effects of the photometric redshift 
systematics on constraints on some parameters in Fig. 8 
and Table VII. 

In Fig. 8, we show the total effects of the scatter and 
bias on the parameter determination, as well as the indi- 
vidual effects. The scatter gives a stronger impact than 
the bias, and including both of the two systematics dra- 
matically enlarges the constraint contours, except for Ji, 
and /nl- Perhaps, for f^, the constraint comes mainly 
from CMB and CMB lensing, so the effects of system- 
atics in the LSS observables do not contribute to the 
constraint. For the constraint on /nl, the effects of sys- 
tematics have little degeneracies with the effect of /nl 
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FIG. 7: Projected la"(68%) confidence constraints in some 
parameter spaces for three fiducial values of the mass thresh- 
old Mobs[/i"^M0] = 10"■^10"■^10'2■^ We show the con- 
tours expected from combining HSC, Planck and ACTPol, 
and these results include all the information from the auto 
and cross correlations. Only in this figure, we plot the devi- 
ation from the fiducial values A log Mobs instead of log Afobs 
for the purpose of illustrating the different fiducial models on 
the same planes. 



because we can distinguish these effects through the ob- 
servation of the power spectra at large-scale region (low-^ 
region) . 

In Table VII, wc show Icr errors for some parameters 
in the case with HSC and Planck experiments, with or 
without ACTPol, sorted by the cross correlation spec- 
trum used in the analysis. The relative amount of con- 
tributions from each cross correlation to the constraints is 
similar to the case without systematics. However, we en- 
joy the benefits of including the extra cross correlations 
more in the case with systematics than without. For 
example, information from the cross correlation brings 
38.7% and 37.7% improvements of the constraints on Oa 
and w in the case with systematics, respectively, while 
bringing 32.6% and 31.3% improvements in the case with- 
out systematics. 



FIG. 8: Projected la"(68%) confidence constraints in some 
parameter spaces without photometric redshift systematics, 
including the effects of scatter and bias individually, and in- 
cluding the both of the effects simultaneously. The contours 
are obtained from HSC, Planck and ACTPol, and the results 
include the information from all auto and cross correlations. 



so surveys observing high-redshift and wide area are pre- 
dictably effective for this purpose. Here, we show the 
constraints expected from other survey projects besides 
HSC, such as DBS [27] and LSST [28] in Fig. 9. 

Comparing the HSC survey with DES, the DES survey 
is more suitable for the constraints on the parameters 
related to dark energy, such as f^A and than HSC. 
For the constraint on the primordial non-Gaussianity, 
however, HSC is more suitable than DES. This is be- 
cause HSC will observe galaxies in higher redshift than 
DES and it allows us to follow the redshift evolution 
of the scale-dependent bias due to the primordial non- 
Gaussianity, even though the survey area of DES is larger 
than HSC. On the other hand, the constraints from LSST 
are tremendous for most of the parameter determina- 
tions. 



VIII. SUMMARY 



E. Comparison with other various future surveys 

Focusing on the constraints on the primordial non- 
Gaussianity from the galaxy power spectrum or two- 
point correlation function using the scale-dependent bias, 
the result significantly depends on the survey strategies. 
The signature of the primordial non-Gaussianity grows 
with redshift and appears prominently on the large-scale. 



In this work, we estimated errors in the determina- 
tion of cosmological parameters for some future surveys 
with the Fisher matrix method, newly including galaxy- 
galaxy lensing (g7) and CMB lensing-galaxy lensing (V'7) 
cross correlations. In general, the extra information from 
cross correlations allows us to estimate the cosmological 
parameters more precisely. 

As for the constraint on the /NL-logMobs plane, 
galaxy-galaxy lensing (g7) cross correlation improves the 
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FIG. 9: Projected 1(7(68%) confidence constraints in some 
parameter spaces for /i/ 7^ model witfi tfie photometric red- 
sfiift scatter. Tiie contours are for DES (red dasfied), HSC 
(blue dot-dashed) or LSST (green), with Planck and ACTPol 
combined. The results include all the information from auto 
and cross correlations. 



constraint greatly. On the other hand, galaxy-CMB Icns- 
ing (gip) cross correlation can also improve the constraint 
if the ACTPol experiment is included. 

Although galaxy-galaxy lensing (g7) cross correlation 
has a great contribution for improving the constraint 
on /nl, we have to pay attention to the systematics of 
galaxy lensing observation. The photo-z error diffuses 
the observed galaxy distribution of each redshift bin, and 
changes the overall amplitude of the power spectrum. 
This behavior is similar to the change of the param- 
eter logA'/obs- On the other hand, positive (negative) 
/nl enhances (decreases) the amplitude only on large- 
scales. Therefore comparing these two parameters, the 
constraint on log Mobs is more affected by photo-z error 
than that of /nl- 

Here, we accounted only for the photo-z error for sys- 
tematics, but it is known that there are other system- 
atics for the galaxy lensing observables. One of them is 
a complication in the measurement of the shear due to 
the incomplete removal of the effects of the point spread 
function (PSF) [65]. So the contribution of galaxy-galaxy 
lensing (g7) cross correlation to the constraint on /nl 
may be diminished if we consider a more realistic con- 
dition including other systematics for galaxy lensing ob- 
servables. In such a case, the impact of galaxy-CMB 
lensing (g?/)) cross correlation probably shines out more. 

In general, the choice of the fiducial cosmological model 



affects the parameter forecasts. In this paper, we consid- 
ered the differences of the fiducial model in the neutrino 
mass and the photometric redshift scatter. As far as the 
constraint on /nl is concerned, we found that the selec- 
tion of the fiducial model makes an impact on the result. 

It should be emphasized that combining the satellite 
CMB experiment (Planck) with the ground based one 
(ACTPol) can greatly improve the constraints on the cos- 
mological parameters. This is because the information of 
CMB lensing can be extracted more efficiently by com- 
bining Planck with ACTPol than by Planck only. For 
the constraint on /nl, the marginalized error can be im- 
proved from A/nl ^ 5.1 to 4.8 by combining ACTPol 
with Planck, because in this case galaxy-CMB lensing 
cross correlation (gV') starts to play an important role in 
the improvement of the statistical error of /nl- 

As for the strategy of the survey, we found that HSC is 
preferable to DES for the constraint on /nl because the 
former can probe higher redshift. We expect A/nl ~ 5.5 
with the DES survey while A/nl ^ 4.8 with HSC. How- 
ever, we never forget the importance of wide field surveys 
to see the signature of the primordial non-Gaussianity 
appearing on large-scale. 

As for the benefits from the cross correlations, we 
found that the cross correlation between CMB lensing 
and galaxy distribution improves the constraints on /nl 
from 5.4 to 4.8 (All.1%), and from 1.1 to 1.0 (A8.3%) in 
the case with LSST. The relative improvement of /nl by 
including cross correlations is less distinct in LSST, but 
the cross correlation will be still significant in the future 
surveys. 
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